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Abstract 

We address the problem of including Coulomb distortion effects in inclusive quasi- 
elastic (e, e') reactions using the eikonal approximation. Our results indicate that 
Coulomb corrections may become large for heavy nuclei for certain kinematical 
regions. The issues of our model are presented in detail and the results are compared 
to calculations of the Ohio group, where Dirac wave functions were used both for 
electrons and nucleons. Our results are in good agreement with those obtained by 
exact calculations. 
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1 Introduction 



Nucleon knockout by electron scattering provides a powerful probe of the dy- 
namics of nucleons in the nuclear medium. The transparency of the nucleus 
with respect to the electromagnetic probe makes it possible to study the entire 
nuclear volume. For light nuclei, the weakness of the electromagnetic interac- 
tion allows for a separation of the soft Coulomb distortion of the electron 
scattering process from the hard scattering event in which, to a very good 
approximation, a single virtual photon transfers energy and momentum to 
the nuclear constituents. Since the kinematic conditions of electron scatter- 
ing can be varied easily, different aspects of the reaction mechanism can be 
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tested. Under conditions in which a single nucleon receives most of the energy 
and momentum transfer, the quasifree electron-nucleon scattering process is 
emphasized. 

Nucleon knockout from heavier nuclei is used to measure the single-nucleon 
spectral function in complex nuclei. At low excitation energy, exclusive mea- 
surements to discrete states of the residual nucleus provide information on 
quasiparticle properties, such as binding energies, spectroscopic factors, spread- 
ing widths and momentum distributions. The specific models used to extract 
the structure information can be tested by comparing parallel with nonparallel 
kinematics or by varying the ejectile energy. The spatial localization of specific 
orbitals also provides some sensitivity to possible density-dependent modifi- 
cations of the electromagnetic properties of bound nucleons. The azimuthal 
dependence of the cross section and the recoil polarization of the ejectile can 
provide detailed tests of the reaction mechanism, which may be useful in delin- 
eating the role of two-body currents or in testing off-shell models of the current 
operator. Measurements at larger missing energy can provide information on 
the deep-hole spectral function or on multinucleon currents. Measurements at 
large missing momentum are sensitive to short-range and tensor correlations 
between nucleon pairs. 

For heavier nuclei Coulomb corrections (CC) may become very large and af- 
fect the measured cross sections; this needs to be accounted for, if one aims 
at a quantitative interpretation of data. Here, we concentrate on the inclusive 
quasielastic scattering process (e, e'), where only the scattered electron is ob- 
served. We will model this process as a knockout reaction where the nucleons 
are hit by the virtual photon emitted by the scattered electron. 

Inclusive scattering provides information on a number of interesting nuclear 
properties: 

-The width of the quasielastic peak allows a dynamical measurement of the 

nuclear Fermi momentum [1]. 
-The tail of the quasielastic peak at low energy loss and large momentum 

transfer gives information on high-momentum components in nuclear wave 

functions [2]. 

-The integral strength of quasielastic scattering, when compared to sum rules, 
tells us about the reaction mechanism and eventual modifications of nucleon 
form factors in the nuclear medium [3]. 

-The scaling properties of the quasielastic response allows to study the reac- 
tion mechanism [4]. 

-Extrapolation of the quasielastic response to A = oo provides us with a very 
valuable observable for infinite nuclear matter [5]. 

For heavier nuclei, these questions obviously can only be addressed once the 
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Coulomb distortion of the electron waves is properly dealt with. 

In this paper, we present an approximate treatment of electron CC for in- 
clusive quasielastic (e, e') reactions modeled as a nucleon knockout process. 
In the plane- wave Born approximation (PWBA), the electrons are described 
as plane Dirac waves, which is a poor approximation for heavy nuclei with 
strong Coulomb fields. In a better approach, called eikonal distorted wave 
Born approximation (eDWBA), we use electron waves which are distorted by 
an additional phase and a change in the amphtude. This phase shift and the 
modification of the amplitude account for the enhanced momentum and a fo- 
cusing effect which occurs when the electron approaches the strongly attractive 
nucleus. 

Calculations with exact Dirac wave functions have been performed by Kim et 
al. [6] in the Ohio group and Udias et al. [7,8]; the present eikonal approxi- 
mation has the advantage that it is relatively simple to implement and that it 
avoids large computational costs. The eikonal method and its higher approx- 
imations, which can be obtained from an iterative procedure, is expected to 
lead to an asymptotic, rather than a convergent expansion [9]. Therefore, its 
use for the calculation of exclusive cross sections may be problematic accord- 
ing to Giusti and Pacati [10,11], but the good agreement for the inclusive cross 
section with the exact calculations by Kim et al. seems to justify the method 
in this case. Additionally, the eikonal approximation used in [10,11] is not 
equivalent to the approach used in this paper, since we do not calculate the 
electron wave phase shift from an expansion around the center of the nucleus. 
For exclusive cross sections, not only exact wave functions for the electrons 
are needed, but also an accurate description of the proton or neutron wave 
functions. In the inclusive case the fine details of the nuclear structure and fi- 
nal state interaction are not important, and the description of the unobserved 
knocked-out hadron can be rather cursory. 

Various approximate treatments have been proposed in the past for the treat- 
ment of CC [10,12,13,14,15,16,17,18,19], and there is an extensive literature 
on the eikonal approximation [20,21,22,23,24]. In particular it has been shown 
that at lowest order, an expansion of the electron wave function in aZ leads 
to the well known effective momentum approximation (EMA) [9], which is 
explained in Sect. 2.2. The EMA has been used to correct Coulomb distor- 
tions for elastic electron scattering data or inelastic data to low lying states, 
although exact descriptions are readily available. For quasielastic (e, e') scat- 
tering a comparison of EMA calculations with numerical results from the 
'exact' DWBA calculation [25,26] indicates a failure of the EMA [27]. 

We point out that the model of the nuclear structure used in this paper is rel- 
atively simple, since our focus is mainly on the electronic part of the problem, 
i.e. the influence of the Coulomb held on the electron wave function. Still, our 
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simple model leads to satisfactory values for inclusive cross sections which are 
sufficient for our purpose. We also find that the EMA fails to account for CC 
for quasielastic (e, e'). This observation is discussed in detail in Sect. 3. 



2 Models and approximations 

2.1 Effective momentum approximation 

The accurate description of the electrostatic potential of the nucleus is an 
important ingredient for the calculation of CC. In a simplified classical picture, 
the electron hits mostly the outer region of the nucleus. This implies that for 
heavy nuclei like lead with a radius of approximately 6 fm, the kinetic electron 
energy is increased in this region by about 20 MeV, which is a non negligible 
modification when the electron energy lies in the range of a few hundreds 
of MeV. In the most naive picture, the effective momentum approximation, 
one typically assumes that the nucleus is a homogeneously charged sphere 
with equivalent radius R^- For highly relativistic particles in a potential V 

— * 

with asymptotic momentum k, one may neglect the mass of the particle d^' — 
y|, |A; I ^ m), such that the energy-momentum relation reduces to = c = 1) 

{E -Vf = k'^ + ^ E -V = k,k=\k\. (1) 

Then the momentum shift A/c of a highly relativistic electron in the region 
where it interacts with the nucleus follows from the potential energy of the 
electron inside the nucleus (C = 1 at the surface or C = 3/2 in the center) 

M^C^, C = 1...3/2, kff^kij + M. (2) 

The standard EMA uses effective momenta corresponding to the central value 
of the Coulomb potential (C=3/2). The EMA cross section of the considered 
process is then calculated by using the plane wave approach, but with the 
electron momenta replaced by their corresponding effective values. In this 
way one accounts for the fact that the electron wave length is reduced in 
the relevant nuclear region. But due to the attractive Coulomb potential, the 
modulus of the initial (final) electron wave is also enhanced by a factor of 
~ Fj {Ff) inside the nuclear volume. The cross section is therefore multiplied 
additionally by a factor Ff = {kf^/k,i)'^ which accounts for the focusing of 
the incoming electron wave in the nuclear center. The cross section is not 
multiplied by Fj = {k'f/kf)^, because this factor is already contained in the 
artificially enhanced phase space factor of the outgoing electron. 
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In order to be more explicit, we mention that in the plane wave Born approx- 
imation, the cross section for inclusive quasielastic electron scattering can be 
written by the help of the total response function Stotai as 

^^^f^ - ^Mo« X S^otaliWl^, ©e), (3) 

ailfaef 



where the Mott cross section is given by 

aMott = ^a'cos\Qe/2)e}/q^^. (4) 



Here, a is the fine-structure constant, Og is the electron scattering angle, 
q = ki — kf is the momentum transfer given by the initial and final electron 
momentum, a; = ej — e/ is the energy transfer given by the initial and final 
electron energy, and the four-momentum transfer squared is given by — 
oj^ — q^. 

The Mott cross section remains unchanged when it gets multiplied by the EMA 
focusing factors and the momentum transfer gr^ is replaced by its corresponding 
effective value, i.e. Fj^F|/g^gg = Therefore, the EMA cross section can 
also be obtained from (3) by leaving the Mott cross section unchanged and by 
replacing Stotai^q\,uj, Ge) by the effective value Stotai^kf -kf\,LO,Q)^. 



2.2 Eikonal approximation 



We can take the local change in the momentum of the incoming particle with 
momentum — kiki into account approximately by modifying the plane wave 
describing the initial state of the particle by the eikonal phase xi (r) 



where 

z 

Xi(0 = - / V{r + hs)ds^- I V(x,y,z')dz' (6) 



if we set ki — k\ez- As desired, the z-component of the momentum then 
becomes 



= -id^e^^^^^'^^ = {ki - v)e'^'^^^'^\ (7) 
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The final state wave function is constructed analogously 



.ikff-ix2ir) 



(8) 



where 



oo 



X2(r) ^- I V{r + kfs')ds' 



(9) 







For the sake of simplicity, we consider spinless electrons in Sects. 2.1-2.4. In 
our actual calculations spin is included. The spatial part of the free electron 
current (which interacts via photon exchange with the particles inside the 
nucleus) 



where e is the charge of the electron. The spatial part of the electron current 
now contains the additional eikonal phase, and the prefactor contains gradi- 
ent terms of the eikonal phase which represent essentially the change of the 
electron momentum due to the attraction of the electron by the nucleus. 

So far we have only considered the modification of the phase of the wave 
function, and for many applications this is a sufficient approximation. It has 
been applied to elastic high energy scattering of Dirac particles in an early 
paper by Baker [28]. However, the method leads, e.g. for quasielastic scattering 
of electrons on lead with initial energy = 300 MeV and energy transfer 
uj = 100 MeV, to errors up to 50% in the cross sections. The reason is that 
also the amplitude of the incoming and outgoing particle wave functions is 
changed due to the Coulomb attraction, as mentioned above. This fact can 
be related to the classical observation that an ensemble of negatively charged 
test particles approaching a nucleus is focused due to its attractive potential. 
Below, we give a simple classical treatment of this fact which illustrates the 
basic properties of the phenomenon. A semiclassical discussion of the focusing 
factor can also be found in [9]. 



is replaced by 



J EIK 




-ikff+ix2 (r) \ ik^f+ixi (r) 



e{h + kf + Vxi - Vx2)e'(^*-^^)''+'(^i+^^) , 
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Nucleus 

Fig. 1. Electrons incident on an nucleus with impact parameter 60. 
2. 2. 1 The focusing factor 

The focusing factor can be derived approximately from a classical toy model 
according to Fig. (1). We consider the trajectories of an ensemble of highly 
relativistic particles approaching the nucleus in the center of the coordinate 
system in z-direction with an impact parameter between 60 and 60 + dh^. The 
longitudinal velocity of the particles can be taken as the speed of light, since 
the particles are highly relativistic and the change of the velocity in transverse 
direction causes only a second order effect to the longitudinal component. 
Therefore we have z{t) = t, r{t) ■ 



'bg + t'^, keeping in mind that the impact 
parameter can be considered as constant at this stage of our approximation. 
The density of the particles at z is increased by a focusing factor / which 
is given by the ratio of the area of two annuli with radii bo, bo + dbo and 
b{bo,z),b{bo + dbo,z): 



r-i,, . db(bo, z) b(bo, z) 



(12) 



We will calculate now b{bo,z) for the screened potential 

aZ 



(13) 



since in this case a simple study of the problem and its most important prop- 
erties is possible. 



The force Ft acting on the particle in transverse direction is given by 
b dV{r) aZbo 



Ft — — 



r dr 



+ i?2 



■3 ' 



(14) 



7 



whereas the "transverse mass" is given by the energy E (c— 1). Therefore we 
obtain for the transverse acceleration 



aZ bo , . 

vt^--^— s , (15) 

^ ^f' + bl + R^ 



and 

aZ bn t + Jt^ + bl + R^ 

vxit) = 



E bl + R^ + bl + R^ 



aZ bo t+Vr^ + R^ 

E bl + R^ r^ + R^ ' ^ ^ 



From VT{t — > oo) = we obtain for a pure Coulomb field {R -C 60) 

the well-known transverse momentum transfer 



Akr = -J—. (17) 



Furthermore, we obtain from (16) 
hf \ K aZboiz + V^TW) 

Kz) = bo-^ 6fT^^ ' ^is) 



The focusing factor is then given by (/3 = aZ/E) 



' Z + ^Jr^ + R^ _ 2bl{z + ^r^ + R'^) bl ^ 

^ bl + R^ {bl + i?2)2 + (^,2 + R2^^/:fT^, 



Since our calculation is not exact, we keep only the relevant first order in /3 
in (19). For the focusing factor in the center of the nucleus follows 

m = (i-^r^(i-^f~(K/k.f, (20) 



i.e. one obtains the focusing factor used in the effective momentum approxi- 
mation, where the increased density near the nucleus is taken into account by 
multiplying the wave function by a suitable factor ~ /^/^. In (20), ki = E/c 
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denotes the asymptotic momentum of the highly relativistic incident particles, 
whereas A;- is the momentum of a particle with impact parameter 60 = in the 
center of the nucleus, which is given by /c^ = fcj + aZ/{Rc). c is the speed of 
light. The amplitude of the wave function is then correctly normalized in the 
region of the nucleus, where the nucleon knockout process is taking place. The 
discussion given above is only classical, but it is sufficient to show the most 
important features of the focusing effect. E.g., from (19) is it obvious that the 
amplitude of the wave function continues to increase on the rear side of the 
nucleus due to the deflection of the incoming particle wave. 

Knoll [14] derived the focusing effect from a high energy partial wave expan- 
sion, following previous results given by Lenz and Rosenfelder [13,18]. For the 
incoming particle wave expanded around the center of the nucleus he obtained: 

{1 + air^ - 2a2k'f+ iair^k'f+ ia2[{k' x r)^ + a{k' x r)]}uk, (21) 



where k' is parallel to k, and where a acts on the spinor Uk to describe spin 
dependent effects, which are very small in our cases of interest. An analogous 
equation holds for the distortion of the outgoing wave. The parameters 01,2 
depend on the shape of the potential. For a homogeneously charged sphere 
with radius Rg they are given by 

aZ 3q:Z 
"^ = "6F^' ""'^'Ik^: ^^^^ 



The increase of the amplitude of the wave while passing through the nucleus is 
given by the — 2a2/c'r-term. Taking our classical result for a screened potential 
leads to 02 ~ aZ/k''^R^, a result which cannot be compared directly to (22) 
due to the different shape of the potentials. Due to the fact that the amplitude 
of the wave function is smaller in the upstream side of the nucleus but larger 
by a similar amount on the downstream side, the influence of the — 2a2A;'r-term 
on the cross section in general is not very large. Minor effects can be observed 
for large scattering angles. Also the spin related term is of minor importance 
for highly relativistic electron energies. 

The Oir^-term accounts for the decrease of the focusing also in transverse 
direction. For the cases of interest in this paper, it leads to 1-2 percent effects 
in the cross sections. Imaginary terms like ia2{k' x r)^ describe the deformation 
of the wave front near the center of the nucleus. They could be obtained 
correspondingly by an expansion of the eikonal phase in that region, and 5+ 
is the eikonal phase in the center of the nucleus. 
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In the present work, the cikonal phase is obtained directly from an analytic 
expression as described in Sect. 2.2, whereas the focusing is calculated using 
Knoll's results given above for a homogeneously charged sphere. 



2.3 Electrostatic potential of the nucleus 



For our present eikonal calculations, we use a potential energy of the electron 
of the form 



i-^rnr) - -^^^2 - ^J^f^. > (23) 



2 + 24 R'^R'r^ 



which goes over into a Coulomb potential for r — oo, and which is a good 
approximation for the potential generated by the Woods-Saxon-like charge 
distribution of a nucleus (see Fig. 2). R' can be used as an additional fit 
parameter. A good choice is R' = 0.5174it!e. Furthermore, expression (23) 
has the advantage that it is possible to derive an analytic expression for the 
eikonal phase, which is convenient when the eikonal phase has to be calculated 
numerically in a computer program. 

The eikonal phase turns out to be divergent for a Coulomb-like potential, 
but it is possible to regularize the eikonal phase by subtracting a screening 
potential ~ (r^ + a^)^^/^ with a ^ R from (23), such that the potential falls 
off like for large r. The divergence can then be absorbed in a constant 
divergent phase ~ log(a) without physical significance, when the limit a — > oo 
is taken. It is quite instructive to calculate the eikonal phase for the simple 
screened potential [29] 

-rw \ (^Z . aZ aZ 

Vi(r) = , , K"(r) = , + , . (24) 



One obtains for a particle incident parallel to the z-axis for impact parameter 
b (r'2 = b^ + z'^) 



} 1 1 

xl^oiZ / dz' ( , - ) = 



(^ + \/r2 + a2)(62 + i?2)' ^ ' 
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where r — b + z , and therefore for the regularized eikonal phase 



^1 = Jl^(^i - «^log(«)) = Q^^log( y^^, )■ (26) 



which is defined up to a constant phase. Taking the gradient of Xi in transverse 
direction 

dxi ^ aZb z + Vr^ + R^ 



gives for the transverse momentum transfer the same result as the classical 
expression (16). This illustrates the fact that the eikonal approximation also 
accounts for the transverse modification of the particle momentum. The charge 
density 

p(r) = -^d'^(rV{r)) (28) 
er 



corresponding to the potential given by Eq. (23) satisfies 

{p)^eZ, (29) 

{r'p) = lR'eZ, (30) 



i.e. we can indeed identify R^ with the equivalent radius of a homogeneously 
charged sphere Rl which is given approximately by 

Re = 1.128A^/^ fm + 2.24A-^/^ fm (31) 



for nuclei with A > 20. Rg can be related to the rms radius Rm by 

Rl - ■ (32) 



For a simple potential ~ (r^ + R^) ^1"^ ^ the rms radius does not exist, since 
the corresponding charge distribution does not fall off fast enough. 

The expressions necessary for the calculation of the eikonal phase for the 
potential (23) are given in the appendix. 
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Fig. 2. Comparison of different model potentials: Potential V used in our calcula- 
tions, potential generated by a homogeneously charged sphere with radius Rg 
and potential generated by a corresponding Woods-Saxon charge distribution. 



2.4 Scattering cross section 



For the sake of notational simplicity and in order to give a transparent in- 
troduction to the method used for the calculation of the CC, we will neglect 
the dependence of the interactions on the spin and internal structure of the 
particles. Therefore, the particle wave functions are replaced by scalar fields in 
the following discussion. We call the positively charged particle the 'proton', 
whereas the negatively charged particle is the 'electron'. The simplified scalar 
expression for the electron current can then be replaced by the correct Dirac 
current (43) in a straightforward way. The nucleon current which we used for 
our calculations is discussed below. 

The lowest order transfer current of an electron with initial/final state wave 
function (f)'^:^ and a scalar (point-like) proton with initial/final state wave func- 
tion 4>'!^^ is given by 

= +ze{<Pi^dyi - <PldyfJ) , (33) 
j^^-ie{<f>pd^<f>i-<l>;d^<f>f;). (34) 



For the scattering cross section of an electron with initial and final momentum 
ki = {ei,ki), /cj = {ef,kf) off a proton with momenta = {Ei,pi), p'j — 
{Ef,pf) we obtain in Born approximation the second order contribution in 
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the coupling constant e to the differential cross section (e^ = Ana) : 

X {2nyS^^\ki+pi-kf-pf)\M\^, (35) 



where the matrix element M is given by the product of the currents and the 
photon propagator 

M = e^(A±MB(El±Pfr. (36) 
[ki — kf) 



The corresponding expression for the physical situation where electrons and 
nucleons are treated as particles with spin can be found, e.g., in [8]. 

Since we consider proton knockout by highly rclativistic electrons from a nu- 
cleus at rest, the velocity in the flux term \v'[ — vf\ will be set to the speed of 
hght. 

If the initial proton is in a bound state with wave function 

</.;(f, t) = (27r)-3/2 J d^q^^^y^^-m^ (37) 

the transition probability is replaced by {pf = {E^, q)) 
{2nY5^^\h+pi-kf -pf)\M\^ ^ 
2T:5{ei + E^ - Sf - Ef) 

{27r)^S{ei + Ei-ef-Ef) 

~^'^'-iP{kf+p-'f-ki)\ (38) 



X 



^(ki + kf)^(2pf + kf-ki)''^.^^ . .^,2 



{ki - kfY 



In the distorted wave Born approximation, the presence of the Coulomb field 

is taken into account by using exact wave functions of the particles in the 
Coulomb field. In the eDWBA, these wave functions are replaced by their 
eikonal approximation. Therefore, the electron current is modified by the 
eikonal phase and due to the focusing of the electron wave. For electrons 
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with potential energy V{r) in the central Coulomb field of the nucleus we 
must replace 



by 



where 

^r = (e,-V^(r),^, + Vxi), (41) 

~^ = [ef-V{r),kf-Vx2) (42) 

correspond to the modified electron current in the eikonal approximation, and 
f{r) denotes the focusing factor resulting from the incoming and outgoing par- 
ticle wave function. The gradient and potential terms in the electron current 
are an artefact of our treatment of electrons as scalar particles. They appear 
in a similar fashion in the Gordon form of the electron current. The standard 
Dirac form of the current 

f = (43) 



can be split into a convective current and a spin current after some algebra 
by making use of the Dirac equation. One obtains 

If — — ^1 P — P — 

f = —\^d^^ - {d^^)^ + — d^\^a^''^\ . (44) 



The convective and the spin current are separately conserved and gauge in- 
variant. But it is not advisable to use the Gordon form (44) in conjunction 
with the inexact eikonal approximation of the electron wave function, because 
the electron current is then described less accurately than by the standard 
expression (43). 

It is now no longer possible to evaluate the coordinate space integral in a 
trivial way. As a first approximation we assume that the initial momentum 
("missing momentum") of the bound proton is given approximately by the 
external momenta, i.e. 

q = kf +Pf -ki- Aq, Ag ~ 0. (45) 
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Then we can expand 




(46) 



where = {ki — kfY 
then given by the term 



^. The zeroth order contribution to (40) is 




J (fr {ki + kf)^{pi +pfY i^ir) f{r)e' 



,i{ki-kf-pf)r+ix{r) 



(47) 



The A g-terms stemming from the expansion above and from the expression 
for the proton current act then as a gradients on the eikonal phase in real 
space according to 



such that we obtain from (46) an expansion of the integral (40) which contains 
higher derivatives of the eikonal phase. For the actual calculations, it proved 
sufficient to include all terms containing first and second order derivatives. 
This way, the matrix clement (40) can be reduced to three-dimensional inte- 
grals which are numerically tractable. This is due to the fact that the integrals 
extend only over a finite volume around the nucleus. For the case presented 
in Fig. 4, second order derivative terms contribute about one percent to the 
total cross section. 

2.5 Nucleon-nucleus interaction 

As a first step, we adopted also an eikonal approximation for the proton wave 
function, in order to improve the poor approximation of describing the outgo- 
ing proton by a plane wave. The corresponding treatment can be carried out 
along the same lines as discussed above. In order to be as realistic as possible, 
we use an energy-dependent volume-central part of an optical model potential 
as given in a recent work [30] . The potential is given as the sum of a Woods- 
Saxon potential and the Coulomb potential of the nucleus for protons. The 



/ 



d^q J d^r AqiP{q)e'^'''-^^-^f-Pf^^+'^^^ = 
(27r)^/'i J dhi){r)e'^^'-^f-^f^''Ve'''^'^, 



(48) 



15 



depth of the Woods-Saxon part of the potential depends on the energy E of 
the proton (in MeV) by 

Vws{E) = vi[l - V2{E - Ef) + vs{E - E^f - Vi{E - Eff] , (49) 
where 

vi = 67.2, V2 = 7.9 • 10"^ V3 = 2.0 ■ 10"^ V4 = 7 ■ 10"^ Ef = -5.9. (50) 

For neutrons we have 

vi = 50.6, V2 = 6.9 • 10"^ ^3 = 1.5 ■ 10"^ V4 = 7 ■ 10~^, Ef = -5.65. (51) 

It has been noted in [31] that when comparing the equivalent central potentials 
generated by the sum of the scalar and vector part of relativistic potentials, 
the real part of the optical potential is 10 MeV deeper for ^°Ca. A comparison 
of a standard nonrelativistic central potential and a corresponding potential 
obtained from a Dirac equation based procedure given in [7] shows good agree- 
ment for 208pjj (^^Y^Q difference is of the order of only 2 MeV), as well as a good 
agreement with the optical model potential given above. 

The imaginary part of the optical potential was not taken into account in our 
calculations. The imaginary part is intended to describe the loss of flux in 

proton (neutron) elastic scattering. In inclusive processes, only the electrons 
arc observed, and one does not need to take into account whether or not the 
ejected nucleons got "lost", i.e. initiated some subsequent nuclear reaction. 

The same wave functions for the outgoing nucleon are used for the PWBA 
and eDWBA calculation, since we are only interested in the electronic part 
of CC. The eikonal approximation, when also used for the calculation of the 
wave function of the outgoing nucleons, leads to a considerable additional 
calculational effort. We found that an effective momentum approximation for 
the proton calculated from the energy-dependent optical model potential leads 
to equally good results, and reduces the calculational effort considerably. It 
is obvious that the plane waves used as final states arc not orthogonal to the 
initial bound states, but the results obtained for the cross sections are still 
quite satisfactory for our purposes despite of the fact that we use a single 
particle and a plane wave approximation. We studied also the inclusion of a 
focusing factor for the proton wave function in order to improve upon the plane 
wave approximation. Although such a procedure might be used to improve the 
theoretical cross sections artificially, it is only an ad hoc procedure and was 
therefore not used for the calculations presented in this paper. 
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Fig. 3. Charge density of the Pb nucleus derived from the harmonic oscillator model 
and a Wood-Saxon potential fit. 

In a first attempt to calculate CC for a heavy nucleus like ^°^Pb, we used 
the harmonic oscillator shell model wave functions for the bound protons. 
These wave functions are easily accessible for numerical calculations due to 
their simple analytic form. Harmonic oscillator wave functions are satisfactory 
approximations for nucleon wave functions especially when they are scaled 
according to the root mean square (rms) radius of the individual shells; due 
to the requirement of orthogonality of initial states, the same scaling must be 
applied to all wave functions. The contribution of partially filled upper shells 
to the inclusive proton knockout cross section was taken into account by a 
weighting factor for the total shell contribution according to the number of 
occupied states. 

As an alternative we also used wave functions generated from a Woods-Saxon 
potential for the protons. The Woods-Saxon potential, which included an LS 
coupling term, was optimized in such a way that the experimental binding en- 
ergies of the upper proton shells and the rms of the nuclear charge distribution 
were reproduced correctly. A comparison of calculated results showed that the 
ratio (Jpy/ba/ (^cc foi' a kinematic situation as presented in Fig. 4 differs only 
by about one percent for the two models. A comparison of the charge distri- 
bution in a Pb nucleus resulting from the two models is shown in Fig. 3. The 
harmonic oscillator model overestimates the charge density in the center of the 
nucleus, but the impact of this mismatch is reduced due to the small volume 
of the central region. The rms radius of the neutron density distribution was 
taken from a recent study presented in [32]. The binding energies, which also 
enter in the calculation of the wave functions of the outgoing protons and in 
the phase space factors, were taken from [33]. 
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The current of a spinless proton could be modeled by the widely used electric 
Sachs form factor Ge [34]: 

Ge{Q'') = (1 + gV(843 MeV)2)-2 , (52) 



which can be expanded by the same method as the photon propagator (46). 
It is expected that (52) also provides a good description of the form factor of 
a proton inside a nucleus [3]. For the energy range considered in this paper, 
the contribution from magnetic scattering must not be neglected. For the 
calculations, we used the free nucleon current given in [6] and introduced by 
de Forest [35] 

Jii = e*n,/ J^*n,i , (53) 



where 



4^(cc2) = Fa" + ^^'^'^-^2, (54) 



and = pI —pv denotes here the difference between the proton four-momenta 
of the initial and final state, k the anomalous magnetic moments of the cor- 
responding nucleon. In configuration space, the three-momenta have to be 
considered as gradient operators. Fi 2 are the nucleon form factors related in 
the usual way to the electric and magnetic Sachs form factors. 

The relativistic form of the wave function for the bound protons and neu- 
trons was constructed from the nonrelativistic wave function. Whereas the 
large component of the corresponding correctly normalized Dirac spinor was 
obtained directly from the nonrelativistic wave function, the small component 
was derived from the large component in a straightforward way by treating 
the nucleon wave function as a free Dirac plane wave. The momenta of the 
nucleons which were used in the current can be calculated from energy and 
momentum conservation. But in order to improve our model for the current, 
we adopted additionally an effective momentum approximation for the nu- 
cleons, where the influence of the energy-dependent optical potential on the 
ejected nucleons was taken into accoimt. We found that different choices for 
the nucleon current do not affect significantly the ratio of the PWBA and 
eDWBA calculations. This is not the case for the cross section themselves. An 
alternative choice to the cc2 current is the ccl current given by the operator 
[35] 

J^iccl) = (Fi + kF,)^^ - U f,, (55) 
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which was also used in order to check how the ratio of the cross sections 
behaves for different models of the nuclear current. None of the expressions 
for the current (ccl,cc2) is fully satisfactory and both expressions fail to fulfill 
current conservation, but given the fact that we focus mainly on the electronic 
part of the problem, the simple choices given above provide a satisfactory 
description of the proton cTirrcnt. For the cc2 choice, integration by parts 
allows one in a simple way to get rid of the momentum operators acting on 
the nucleon wave functions. 

Some critical remarks which highlight the approximations made within our 
single particle model are in order. Due to the fact that we are working in 
real space, we scaled the nucleon wave functions according to the experimen- 
tal rms radius of the nucleus. This does not automatically imply that the 
momentum distribution of the nucleons is described very accurately. The dif- 
ference between the theoretical Fermi momentum obtained in our model and 
the experimental value is of the order of 10%. Additionally, the choice to use an 
effective momentum approximation for the nucleons is an ad hoc prescription, 
which differs from the naive plane wave approximation for nucleons where the 
influence of an optical potential is neglected. In both cases, one has to accept 
that focusing effects of the nucleon wave functions in the nuclear region are 
missing and furthermore that unitarity is violated to a certain degree. 

Kim et al. [6] presented calculations where plane waves and Dirac wave func- 
tions were used for the electrons, but the nucleon wave functions were cal- 
culated in both cases from a relativistic a — uj model. For the kinematical 
situation presented in Fig. 6 in the next section, their theoretical cross sec- 
tions are slightly below the experimental values, whereas our model leads to 
cross sections slightly above the experimental values. The advantage of our 
model is that we obtain a curve which shows a very similar behavior as the 
experimental curve, although this should not be considered as a virtue of our 
method. We expect that our model provides reliable results only for the ra- 
tio of cross section (were one uses plane waves or 'exact' wave functions for 
the electrons), which can be used for the analysis of experimental data. We 
checked this assumption by varying the optical potential, the binding energies 
and by rescahng the wave functions within reasonable limits. The relatively 
large impact of such variations on the cross section is divided out for the most 
part in the ratio of the cross sections. 

Kim et al. [36] presented also a simplified model where free plane wave func- 
tions for the ejected nucleons and harmonic wave functions for the bound 
nucleons were used. Adopting the same strategy in our case, we obtain nearly 
identical results as those presented in [36] when we scale our wave functions 
according to the experimental value of the nuclear Fermi momentum. 
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The numerical evaluation of (40) was performed by putting the nucleus on 
a three dimensional grid with a side length of 36 fm and using the Simpson 
method as a very simple but efficient integration tool. The number of neces- 
sary grid points is mainly dictated by the wave length of the oscillatory term 
^i{ki-kf-pf) g^j^^ ^g^g ^Yie range of 60^ to 120^ in order to ensure an accuracy 
of 10~^ percent for the values of the integrals. 

It is instructive to have a (classical) look a the approximate size of the different 
effects which lead altogether to the CC. As a specific case we choose the initial 
energy of the electron = 485 MeV, scattering angle 9e — 60° and and energy 
transfer u = 100 MeV. Additionaly, we take the Coulomb potential energy 
V^(0) of an electron in the center of a ^'^^Pb nucleus as —25 MeV. The focusing 
then enters in the cross section as a factor (510/485)^(410/385)^ ~ 1.25, i.e. 
it leads to an increase of the cross section of about 25%. Prom the photon 
propagator term ~ on the contrary we obtain a reduction by a factor 
of 0.80. If we use the potential at the surface of the nucleus —17 MeV, the 
reduction factor is 0.85. This discrepancy illustrates the importance of having 
an accurate description of the electron wave function. Furthermore, the form 
factor of the proton enters the cross section as a factor of the order of ~ 
(1 + QV(843MeV)2)-4, reducing the cross section by 6-9% (for V(0)= -17 
Mev to —25 MeV). Finally, the nontrivial change of the interaction which 
involves the detailed structure of the current via the nucleon wave functions 
gives an effect of several percent. All these effects described above combine to 
a CC which is almost zero for the present example. 



3 Results 

3. 1 Comparison of the results to other approaches 

Fig. 4 shows the ratio of the inclusive cross section for nucleon knockout 
with and without CC for electrons with initial energy of = 485 MeV and 
scattering angle ©e = 60° on ^°^Pb. Neutron knockout contributes by about 
20% to the total nucleon knockout cross section in the considered kinematical 
region. The results obtained by Kim et al. [6] and the eikonal approximation 
are in good agreement, also for different kinematics as e.g. shown in Fig. 5. 
The result for the local effective momentum approximation (LEMA) of Kim 
et al, which is based on parameters that were fitted to the exact calculation, 
is also fairly close. 

The conditions presented in Fig. 5 correspond to a similar momentum transfer, 
but large scattering angle Qg = 143° and incident energy 310 MeV. The initial 
electron energy is smaller, and the effect of CC is more pronounced. CC lead 
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Fig. 4. Comparison of the Coulomb corrections for different approaches. The dotted 
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Fig. 5. c^pwbaI'^cc fo'^ — MeV and 6e = 143°. Dotted curve: (tpwba/c^ema 
as in Fig. 4. 

to a sizeable modification in the longitudinal and transverse responses that 
can be extracted from the data via Rosenbluth plots. The linearity of the 
Rosenbluth plots served also as an independent check for the validity of our 
plane wave calculations. 

The calculations clearly indicate that the effective momentum approximation, 
which leads to admissible results for light nuclei where the CC are relatively 
small, is useless for highly charged nuclei. The curve in Fig. 5 was calculated 
by choosing a Ak = 25 MeV/c according to Eq. (2). Taking smaller values for 
Ak corresponding to the electrostatic potential at the surface of the nucleus 
does not change the situation significantly. 
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Fig. 6. Comparison of theoretical cross sections obtained in PWBA and eikonal 
approximation to experimental Saclay data (e^ = 310 MeV and 6e = 143"). 



It is important to take into account the final state interaction of the proton 
by an energy-dependent optical potential. For protons with energies above 
100 MeV, the potential is becoming increasingly repulsive. Calculations for 
ijj > 180 MeV in our single particle shell model are not applicable, since 
pion production is not included in our calculations, and correlation effects arc 
becoming increasingly important for large uj. Results for < 80 MeV may 
also become dubious due to the large distortion of the outgoing proton wave 
function by the final state interaction, but the simple approximation given in 
this work allows for a good estimation of CC in a relatively wide range of 
energies. 

It is a well-known phenomenon that plane wave calculations have a tendency 
to overestimate cross sections in the region of low electron energy loss, whereas 
cross sections are underestimated in the high energy loss region. The reason 
for this is the fact that the use of plane waves for the ejected nucleons is clearly 
a rough approximation, and the generally used single particle model does not 
contain contributions arising from correlations and meson-exchange currents. 
Inelastic scattering from the nucleons, i.e. excitation of the delta resonance, is 
also absent. 

Nevertheless, the agreement between experimental data taken at Saclay [37] 
and our eikonal corrected calculations is quite satisfactory, as displayed in Fig. 
6 in order to give a typical example. Like in the case ©g — 60", one observes 
that the CC become very large for high energy transfer uj ^ 180 MeV. This 
result should be handled with care, since the calculated cross sections become 
small in this energy region. The electroproduction of the delta resonance will 
have to be included to more reliably estimate the CC at very large energy loss. 
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Overall, we find that the eikonal approach gives results that are very close to 
the results from exact calculations. This approach can thus be used in practi- 
cal applications where one greatly benefits from the much lower calculational 
effort involved in the eikonal description. At the same time, our calculation 
shows again that the often-used effective momentum approximation is inade- 
quate for a quantitative treatment of Coulomb corrections. 



3.2 Study of the Coulomb distortion by comparing quasielastic electron and 
positron scattering 



The comparison of experimentally determined quasielastic cross sections in- 
duced with electrons and with positrons represent an ideal test for any theoret- 
ical approach to quantitatively describe Coulomb distortions. Here we consider 
a recent experiment by Gueye et al. [38] which was used to give an experi- 
mental proof of the validity of the EMA for inclusive quasielastic electron 
scattering. The experiment by Gueye et al. compared quasielastic scattering 
of electrons on ^°^Pb with quasielastic scattering of positrons at the same 
kinematics. 
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Fig. 7. Scaling function F{y) obtained from cross sections measured by Zghiche et 

al. on ^'^^Pb at an electron scattering angle of 60 degrees and initial electron energy 
E. F{y) is obtained from the cross sections by dividing out the Mott cross section 
and the g-dependence of the nucleon form factors (both smooth functions of incident 
energy). 

With the assumption that the EMA is a valid approach the authors interpo- 
late the experimental e^-data linearly to the effective energy of the e+-data, 
Eg+ — Eg- — 2Vc with the Coulomb potential Vc and the relative normalization 
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as fit parameters. They find for Vc a value close to the surface and a normal- 
ization factor which is compatible with the experimental uncertainties. For the 
interpolation the large body of inclusive quasielastic data on ^osp]^ measured 
by Zghiche et al. [37] is used. The authors conclude that a modified EMA with 
a yc{f) — 18.7 MeV close to the surface is a vahd approach to correct cross 
section data. 

However, this conclusion is highly questionable when the quality of the data 
by Zghiche et al. is studied. The peculiar behavior of the data is shown in Fig. 
7, where the scaling function F{ii) (essentially the cross section but with the 
^'-dependence of the nucleon form factor and the kinematical broadening of the 
quasielastic peak removed) is shown (see [3] for detailed definitions) . Although 
the electron energies progress in regular steps between the various data sets, 
the scaling functions F{y) display a curious behavior. The peak values for 
F{ii) for two consecutive energies coincide and jump for the following to a 
significantly lower value. In Fig. 10 in [37], the staircase-hke behavior of the 
data is not visible, since cross sections are plotted, which depend strongly on 
the initial electron energy e^. Additionally, the cross sections for el = 420 
MeV and = 600 MeV are missing there. But the staircase-like behavior has 
a detrimental effect upon the interpolation of the electron data as shown by 
Gueye et al. Also, the e"*" and e~ data of Gueye et al. have been e"^ beam 
omittance corrected and normalized to the e~ data shown in Fig. 7, hereby 
further affected by the staircase behavior. We therefore cannot consider the 
experiment as a proof of the validity of the EMA. 

We performed eikonal and effective momentum calculations for quasielastic 
electron (e) and positron (p) scattering on ^osp]-, fQj. kinematics also used in 
the experiment of Gueye et al, i.e. for initial energy ef = ef = 420 MeV and 
scattering angle Oe = ©p = 60°. The calculations for the positrons were per- 
formed in a strictly analogous way as in the case of electrons. Using an effective 
surface value of 17 MeV for the repulsive (attractive) Coulomb potential in the 
case of the positron (electron) EMA calculation, it turns out that the eikonal 
approximation and the EMA are clearly incompatible. The observed discrep- 
ancy can be explained by the assumption that the EMA underestimates the 
defocusing of the positron waves and the focusing of the electron wave in the 
nuclear volume. 

For the sake of completeness, we also show in Fig. 10 experimental cross sec- 
tions for positrons derived from the total response function given in [38], which 
are in satisfactory agreement with the theoretical values shown in Fig. 9. As 
explained before, the difference between experimental data and the theoretical 
cross sections results from the simple nuclear model employed (lack of high 
momenta in the corresponding spectral functions and the absence of delta 
excitations and meson exchange current contributions). 
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Fig. 8. Theoretical cross sections for inclusive quasielastic electron scattering off 
208p|-, fQj. initial electron energy e| = 420 MeV and scattering angle @e = 60°. 
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Fig. 9. Theoretical cross sections for inclusive quasielastic positron scattering off 
208p|-, fQj. initial positron energy = 420 MeV and scattering angle @p = 60°. 

4 Conclusions 



In the present paper we have investigated the role of the Coulomb distortion in 
quasielastic electron nucleus scattering. A reliable treatment of this distortion 
is needed in particular for a determination of the longitudinal response func- 
tion and for an extrapolation of nuclear responses to infinite nuclear matter. 

We have developed an approximate description -the eikonal approximation- 
that is more transparent and numerically easier to deal with than the exact 
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Fig. 10. Experimental positron cross sections derived from the total response func- 
tion given by Gueye et al, compared to the theoretical values shown in Fig. 9 
obtained from the single particle model. 

treatments (solution of the full Dirac equation) previously employed [6,7,8]. 
At the same time, the eikonal approximation is much more realistic than the 
effective momentum approximation often employed in the absence of results 
from exact calculations. 

We find that the eikonal results for the Coulomb distortion are very close to 
the results of exact calculations. We also find that the effective description of 
the Coulomb distortion gives a rather poor description for high values of the 
nuclear charge Z. 
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Appendix 

The potential energy of the electron according to Eq. (23) can be decomposed 
into three parts V{r) — Vi{r) + V2{r) + Vs{r), where 



V^i(r) 



aZ 



V2{r) 



aZR^ 



(^2 + ^2)1/2' 



2(^2 + ^2)3/2' 
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The regularized eikonal phase Xi generated by Vi is given above by Eq. (26). 
Since ¥2^3 decrease faster than for large r, their corresponding contribu- 
tions to the total eikonal phase for a particle incident parallel to the z-axis 
with impact parameter b 

z 

X2,3(r-) = X2,3(&, ^) = - / <^^' ^2,3(r'), r'^ = + (57) 

—00 



need not be regularized. The integrals in (57) are given by 

2 (62 + i?2) (^2 + ^^2)1/2' y^^) 

257r _(}^^ ^ J_\ (A 7^ 

aZRm'^^ - l^Y^ + 6^ + 26^J arctanj^- j - + ^/2)2 

56S 156S 86^^ 



&'2(r2 + i?'2)3 ^,/4(^2 + ^/2)2 26'6(r2 + i?'2) (^2 + ^/2)3 

262^ 362^ Ab'^z 3z 



6'2(r2 + i?'2)2 5/4(^2 + ^/2) (^2 + ^/2)3 2b'^{r^ + R'^) 

37r564 + 2W + i,/4 
+ T ' ^^^^ 



where = fe^ + z^ and 6' = yW+R!^. The coordinate independent definitions 
of the parameters b, and b' are 2; = {k,ir)/k,i and b'^ = — z"^ , b' = yW+R!^ , 
where hi is the momentum of the incident electron. The eikonal phase of the 
wave function of the outgoing electron with momentum kf, which must be 
added to the eikonal phase of the incoming electron in the matrix element of 
the calculated process, can be obtained directly from the expressions above 
by replacing z — —{kfrj/kf and using again the definitions 6^ = — 2;^, 
6'- V62 + i?'2. 
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